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Abstract 

Analyzing 600 pb _1 of e + e~ collisions at 4170 MeV center-of-mass energy with the CLEO- 
c detector, we measure the branching fraction B(Df — > t + v) = (5.52 ± 0.57 ± 0.21)% using the 
r + — > p + F decay mode. Combining with other CLEO measurements of B{Df — > r + i^) we determine 
the pseudoscalar decay constant fo s = (259.7 ± 7.8 ± 3.4) MeV consistent with the value obtained 
from our Df — > measurement of (257.6 ± 10.3 ± 4.3) MeV. Combining these measurements we 
find a value of fo s = (259.0 ± 6.2 ± 3.0) MeV, that differs from the most accurate prediction based 
on unquenched lattice gauge theory of (241 ± 3) MeV by 2.4 standard deviations. We also present 
the first measurements of B{D+ -» K°tt + tt°) = (1.00 ± 0.18 ± 0.04)%, and B{D+ -» 7r + 7r vr ) = 
(0.65 ± 0.13 ± 0.03)%, and measure a new value for B{Df -» r?p+) = (8.9 ± 0.6 ± 0.5)%. 

PACS numbers: 13.20.Fc, 12.38.Gc, 14.40.Lb 
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I. INTRODUCTION 



The purely leptonic decay of the meson occurs in the Standard Model (SM) via the 
annihilation of the constituent charm quark with the constituent anti-strange quark into a 
virtual W + boson that subsequently materializes as a lepton-antineutrino pair. The SM 
decay rate is given by [l| 



T(Dt - tu) = s m]M Dt (l - Jg- ] \ Vea \' , (1) 

where M D + is the mass, is the mass of the charged final state lepton, Gp is the 
Fermi coupling constant, \V CS \ is a Cabibbo-Kobayashi-Maskawa matrix element, and fp> s 
is the "decay constant," a parameter related to the overlap of the heavy and light quark 
wave-functions at zero spatial separation. 

Measurements of Df —>■ fi + u and Df — > t + v have been made with increasing precision 
recently, and a disagreement has emerged between the theoretical value of fo a computed by 
Follana et al. ^ , and the average of these measurements [3, 0] . It has been pointed out by 
Akeroyd and Mahmoudi jl] that physics beyond the SM could contribute differently to p, + v 
and t + v final states, so increased precision on each of these is important. 

We report here on a new measurement of B(D+ — > t + u) using the r + — > p + V decay mode. 
Previously CLEO has reported on this rate using the r + — > n + V and r + — > e + z/F modes. 



II. EXPERIMENTAL METHOD 
A. Selection of D s Candidates 

The CLEO-c detector ^ is equipped to measure the momenta and directions of charged 
particles, identify them using specific ionization (dE/dx) and Cherenkov light (RICH) [3], 
detect photons and determine their directions and energies. 

In this study we use 600 pb _1 of data produced in e + e~ collisions using the Cornell Elec- 
tron Storage Ring (CESR) and recorded near a center-of-mass energy (E CM ) of 4.170 GeV. 
At this energy the e + e~ annihilation cross-section into D~D* + + D*~D^ is approximately 
1 ill. @. 

In this analysis we fully reconstruct a sample of D~ in nine "tag" modes and then find 
candidate p + — > 7r + 7r° decays in this sample. (Mention of any specific decay implies the use 
of its charge-conjugate as well.) The tagselection is identical to that used in our — > /i + z/ 
paper that can be consulted for details [9J. Briefly, we select candidates using their invariant 
masses. We require that the candidate energies are consistent with those expected of a D s 
or D* in D S D* events. The invariant mass distribution for all tag candidates is shown in 
Fig. [TJa). Then we detect an additional photon candidate from the D* decay, and construct 

MM* 2 = (Ecu - E Ds - £ 7 ) 2 - (p CM - Pd s - P 7 ) 2 , (2) 

where Eqm (pcm) is the center-of-mass energy (momentum), Ed s (pdJ is the energy (mo- 
mentum) of the fully reconstructed D~ tag, and E 1 (p 7 ) is the energy (momentum) of the 
additional photon. In performing this calculation we use a kinematic fit that constrains the 
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decay products of the D~ to the known D s mass and conserves overall momentum and en- 
ergy. All photon candidates in the event are tried, except for those that are decay products 
of the D~ tag candidate. Regardless of whether or not the photon forms a D* with the tag, 
for real D*D S events the missing mass-squared MM* 2 , recoiling against the photon and the 
D~ tag should peak at the mass-squared. 
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FIG. 1: (a) Invariant mass of D~ candidates summed over all decay modes and fit to a two-Gaussian 
signal shape plus a straight line for the background. The vertical dot-dashed lines indicate the ±17.5 
MeV definition of the signal region. No MM* 2 cut has been applied, (b) The MM* 2 distribution 
summed over all modes. The curves are fits to the number of signal events using the Crystal 
Ball function and two 5th order Chebyshev background functions; the dashed curve shows the 
background from fake D~ tags, while the dotted curve in (b) shows the sum of the backgrounds 
from multiple photon combinations and fake D~ tags. The vertical dashed lines show the region 
of events selected for further analysis. 

The MM* 2 distributions for events in the D~ invariant mass signal region (±17.5 MeV 
from the D s mass) are shown in Fig. Hfb). In order to find the number of tags used for 
further analysis we perform a two-dimensional binned maximum liklihood fit of the MM* 2 
distribution and the invariant mass distribution in the interval ±60 MeV from the D s mass 
and 3.50 < MM* 2 < 4.25 GeV 2 . The background has two components, both described by 
5th order Chebyshev polynomials; the first comes from the background under the invariant 
mass peak, defined by the sidebands, and the second is due to multiple photon combinations. 
In both cases we allow the parameters to float. 

We find a total of 43859±936±877 events within the interval 3.782 < MM* 2 < 4.000 GeV 2 
and having an invariant mass within ±17.5 MeV of the D s mass, where the first uncertainty 
is statistical and the second is systematic. 



B. Signal Reconstruction 

We select events with one and only one charged track with opposite sign of charge to the 
tag, that is positively identified as a pion. The event also must contain at least one 7r° — > 77 
candidate with an invariant mass divided by the error on the invariant mass (Pull) < 3; if 
there is more than one such candidate we choose the one with the minimum Pull. Tracks 
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or photons that are used as part of the D~ tag are not considered. Unfortunately, hadron 
tracks can and do interact in the detector material, and deposit additional energy in the 
electromagnetic calorimeter. Thus we do not reject events with more than one 7r° candidate. 
Photon candidates must have an energy deposition in the calorimeter consistent with that 
expected for an electromagnetic shower and deposit more than 30 MeV in the barrel or more 
than 50 MeV in the endcap. In principle for Df — > t + u, t + — > p + V decays all the energy 
should be accounted for in the decay products of the tag and and the p + . We sum up any 
energy in the calorimeter not matched with tag or the p + and call this parameter -E ex tra- 
We also compute the MM 2 as 

MM 2 = (Ecu ~ E Ds -E n - E p f - (p CM - Pd s - p 7 - V P f , (3) 

where E p (p p ) are the energy (momentum) of the candidate p + and all other variables are 
the same as defined in Eq. (j2J). While the MM 2 does not peak at zero, because there are two 
missing neutrinos, it is still a useful variable as two-body Df decay backgrounds will peak; 
e.g., — > p + r] peaks at the rj mass-squared. 

We proceed by defining a 7r + 7r° mass window consistent with the p + mass. Fortunately, 
the mass distribution and the branching ratio for r + — ► p + T> decays are well measured [101 ] . 
We select events within ±250 MeV of the p + mass. This selection is chosen to maximize 
efficiency while still not including too much background. This mass selection is 89.3% efficient 
for t + — ► p + V events. 

The expected MM 2 and E extr& distributions from signal Monte Carlo simulation of D~ 
tag and Df — > t + v events are shown in Fig. [2j Here we have included the above mentioned 
selection window on the 7r + 7r° mass. The MM 2 signal shape is fit to a function that is 
the sum of two bifurcated- Gaussian functions. (A bifurcated-Gaussian shape has different 
widths below and above the mean.) 
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FIG. 2: Monte Carlo generated distributions for signal — > t + v, t + — > p + V (a) MM 2 , and 
(b) -Eextra- The curve in (a) is a fit to the data points with the sum of two bifurcated-Gaussian 
functions. 
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C. Background Expectations 



There are two general sources of background expected arising from either combinatoric 
background in the reconstructed D~ tag sample, or specific decay modes of the Df. The 
former are determined by using sidebands of the candidate D~ invariant mass distribution. 
The latter could arise from modes involving p + decays such as i]p + , but could also come from 
any mode that includes a tt + and a tt°, or a rr + and extra energy that is called a tt°. Previous 
studies have shown that requiring the n + candidate to project to the primary event vertex 
eliminates fake charged tracks as a background source . Our first look at the background 
from Z?+ decays uses Monte Carlo simulation. The background MM 2 and E extTa distributions 
are shown in Fig. [3J The specific background modes are enumerated in Appendix A for three 
different intervals of extra energy, E extia < 0.1 GeV, which we expect is dominated by signal, 
0.1 < -Ecxtra < 0.2 GeV, which we expect has similar amounts of signal and background, 
and 0.8 GeV< -Eextra, where the signal is absent. We separate into these three intervals in 
order to test our understanding of the background. 
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FIG. 3: Distributions for generic Monte Carlo backgrounds from Df decays (solid line) in (a) MM 2 
for Eextra < 0.1 GeV, and (b) E cxtrSL for -0.2 < MM 2 < 0.6 GeV 2 . The expected signal (dashed 
line) and total (dotted line) are also shown. The distributions are normalized to the Monte Carlo 
expectations. 

Three final states cause narrow peaks in MM 2 : (i) K° 7T + 7r° peaks at the K° mass-squared 
and has not been previously measured, (ii) 7r + 7r°7r° peaks at the 7r° mass-squared and is also 
unmeasured, and (iii) r]p + peaks at the r\ mass-squared, and is poorly determined. In order 
to properly treat the background we measure the branching fractions of these modes using 
a double tag technique as described in Appendix B. 



D. MM 2 Resolution 

While the MM 2 resolution of the signal is not an issue because the signal does not form 
a narrow peak, several of the backgrounds do have narrow structures and it is necessary to 
model their shapes properly. In Fig. H] we compare the MM 2 distribution for — > r]p + 
signal from the Monte Carlo simulation with the one found in the data where the rj — > 77 
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decay was detected. We also require that 777r + 7r° mass be between 1.85 and 2.10 GeV, 
in order to have a relatively clean sample but not distort the MM 2 shape. The MM 2 is 
computed using Eq. [3] while ignoring the two photons in the r\ decay. 



1630909-017 




MM 2 (GeV 2 ) MM 2 (GeV 2 ) 



FIG. 4: MM 2 spectrum of r]ir + ir , where r\ — ► 77 was detected, but ignored in the calculation for 
both (a) data and (b) Monte Carlo events. The dotted line in (a) is the signal and the dashed line 
the background determined from the fit to the D~ candidate invariant mass sidebands. The signal 
shapes are the sum of a Crystal Ball and Gaussian functions (see text). 

In both the case of the Monte Carlo simulation and the data we fit the signal with the 
sum of Crystal Ball (CB) [l2| and Gaussian functions. We fix some fit parameters that we 
find from the Monte Carlo simulation including the ratio of the r.m.s. widths (c's) for the 
Gaussian and CB functions, set to a value of 6, and the area of the Gaussian function with 
respect to the CB function to be 20%. In the data fit we also include the background given 
by the D~ invariant mass sidebands. The CB function parameters a and N are taken from 
the Monte Carlo simulation as 1 and 4.5, respectively. We find 

(0.0289 ± 0.0006) GeV 2 for the Monte Carlo CB function, (4) 
(0.0320 ± 0.0020) GeV 2 for the data CB function, 

(0.014 ±0.005) GeV 2 . 

We use the resolution as found in the data above when fitting the data for the r]p + 
component and increase the width of the other narrow components. 



O'MC — 
C"Data = 



III. SIGNAL EXTRACTION 

We proceed by preforming a simultaneous fit to the D~ invariant mass, using a mass 
range within ±70 MeV of the nominal mass and the MM 2 . The procedure is similar to that 
used in Ref. [9J. We first check our procedures by fitting the data MM 2 distribution in the 
-E-extra interval above 0.8 GeV, where we have only background. We include the following 
decay modes as individual probability density functions (PDFs) in the fit: K°tt + tt°, 7r + 7r°7r°, 
07r + , f]7c + , f]p + , ?7 / 7r + , f]'7T + 7T°, U7T + 7i°. All narrow resonance structures are smeared by an 
additional r.m.s. resolution of 0.014 GeV 2 , as determined by our r]p + study (see section |HD]) . 
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The other modes are lumped together into one other PDF. In the likelihood fit we add 
Gaussian constraints on the expected yields based on the known branching ratios and their 
errors. 

The resulting fit to the data for 0.8 GeV < E cxtia is shown in Fig. [51 The fake D~ back- 
ground has been accounted for by simultaneously fitting the sidebands in D~ invariant mass. 
The two-body modes show evident peaks and are well described by the fit, demonstrating 
that our understanding of the backgrounds appears to be adequate. 
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FIG. 5: Fit to the data (points) for 0.8 GeV < E extra . The various components are r]p + (dotted), 
fake D~ (dashed), K ir + ir (long dash), sum of 7r + 7r 7r°, rjir + , W7r + 7r°, cfiir + , r]'ir + , and rj'ir + ir° 
(dash-dot-dot), and other backgrounds (dashed-dot). The solid curve shows the total. 

We next fit the two bins -E ex tra < 0.1 GeV and 0.1 < -Eextra < 0.2 GeV separately. The 
background PDFs are shown in Figs. [6] and [7] for each interval. Note that they now include 
a separate PDF for Df — > r + u, where the r + decays into either ir + T7, or 7i + 7i°7i°V. The fit 
projections for both intervals are shown in Figs. [8]and[9l For the first interval, £"extra < 0.1 
GeV, the background level from all Df decays is about the same as that from fake Dj, 
and both are considerably smaller than the signal. For the second interval the signal and 
background levels are about equal. Table [I] summarizes the signal and background yields. 
(The notation "r+ -> (vr+ + 7r + 7r°7r°)77" denotes the sum of two modes where Dt — > t + u, 
and the r + decays to either 7r + z7 or 7r + 7r°7r°I7.) 

Adding the signal yields in the two ^extra intervals, taking into account the efficiency for 
finding the p + in each interval, and dividing by the number of D~ tags (43859±936±877) 
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TABLE I: Signal and background yields from the fit in two E ext r& intervals. We also list the 
measured or assumed background branching fractions and the r.m.s. error on the fit constraint 
resulting either from the branching fraction error or other considerations. The # MC indicates the 
predicted background number of events for the assumed branching ratio input as the starting point 
of the fit, while # Data gives the number determined by the fit. 



Component 


B(%) 


Constraint 


-Eextra 


< 0.1 GeV 


0.1 < -Eextra 


< 0.2 GeV 






Error (%) # MC 


# Data 


# MC 


# Data 


Signal 








155.2±16.5 




43.7±11.3 


K°7T+7T° 


1.0±0.2 


20 


26.1 


25.2±4.8 


11.0 


10.5±2.1 


VP + 


8.9 ±0.7 


4.2 


7.1 


7.0±0.6 


10.6 


10.5±0.9 


TT+TT TT° 


0.65 ±0.14 


22 


2.8 


2.8 ±0.6 


1.5 


1.6 ±0.3 


T+ -> (7T+ + TT + TT°TT°p 


1.14 ±0.06 


25t 


8.5 


8.4 ±2.1 


12.2 


10.9 ±3.0 




0.576 ± 0.045 


5.4 


1.0 


1.0 ±0.1 


0.48 


0.5 ±0.1 


r]ir + 


1.58 ±0.21 


13.3 


0.9 


0.9 ±0.1 


0.9 


0.9 ±0.1 




4.35 ± 0.35 


8 


1.7 


1.7 ±0.2 


2.8 


2.8 ±0.3 


Xfj,+i> 


5.9 


35* 


3.4 


3.4 ± 1.2 


7.4 


6.6 ±2.6 


Other background 




30* 


11.5 


11.4 ±3.3 


11.8 


10.5 ±3.3 


Fake D~ background 








81.8 ±5.0 




74.8 ± 4.6 



| The error is based on the uncertainties of the resonant substructure that can alter the efficiency. 
t We assign a 35% uncertainty based upon the error on B{Df — > Xe + u). 
* We assign a 30% uncertainty based on the sample size. 
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FIG. 6: Fits to Monte Carlo simulation for the individual background PDFs for -E^tra < 0.1 GeV, 
for the modes (a) i^vr+vr , (b) vr + vr 7r , (c) r?p+, (d) 7]tt + , (e) cpir + , (f) fi+v, (g) D+ -» t + u, 
t + — > (-7T + + 7r + 7T 7r°)!7, (h) sum of the other small modes. 
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FIG. 7: Fits to Monte Carlo simulation for the individual background PDFs for 0.1 < -Eextra < 0.2 
GeV, for the modes (a) K°tt + tt°, (b) 7r+vr 7r°, (c) i] P + , (d) 7/7T+, (e) ^>7r + , (f) p+u, (g) D+ -> t+u, 
t + — > (7r + + 7r + 7r°7r°)!7, (h) sum of the other small modes. 
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FIG. 8: Fit to the data (points) for E cxtTa < 0.1 GeV. The various components are signal (thick 
solid line), r]p + (dotted), fake D~ (dashed), K°ir + ir° (long dash), sum of 7r + 7r 7r , t]tt + , (f)ir + , 
t + — > (7r + + 7r + 7r°7r°)l7, fi + v, and Xfi + u (dash-dot-dot), and other backgrounds (dashed-dot). The 
thinner solid curve shows the total. 
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FIG. 9: Fit to the data (points) for 0.1 < E cxt r& < 0.2 GeV. The various components are signal 
(thick solid line), r/p + (dotted), fake D~ (dashed), K®it + it® (long dash), sum of 7r + 7r°7r°, r/7r + , 07r + , 
t + — > (-7T + + 7r + 7r 7r°)l7, and Xfi + u (dash-dot-dot), and other backgrounds (dashed-dot). The 
thinner solid curve shows the total. 
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we find 

B{D+ -> t + v) = (5.52 ± 0.57 ± 0.21)%, (5) 

where the first error is statistical and the second systematic. We will discuss the systematic 
errors in the next section. In the smallest E ext ra. interval the branching fractions is (5.48 ± 
0.59)%, while in the higher interval it is (5.65±1.47)%. The numbers are consistent. We note 
the data including the background components is well-modeled in all three ^extra intervals, 
confirming our understanding of the background. 

A. Systematic Errors 

The sources of systematic errors in the branching fraction are listed in Table [Til As 
we have let the branching fractions of the background components float in the fit by their 
known errors, there is no additional contribution from this source. The systematic error in 
the background is estimated using two different techniques. First of all, if we remove the 
Gaussian constraint on the sum of the other small mode background fractions we observe a 
1.1% increase in the signal yield. Secondly, if we change the parameters of the background 
shape containing the sum of the other small modes the yield decreases by 0.5%. A separate 
source of error is the efficiency on the detection of background events; if we change the E extra 
efficiency and the ir° efficiency by their errors, and thus change the background yields, we 
observe a combined error of ±1.1%. An additional systematic error could arise from Cabibbo 
suppressed r + — > K + 7r°T7 decays. The measured branching fraction for these decays is 1.6% 
of that of 7T + 7r°I7. The combination of dE/dx and RICH particle identification reduces the 
kaon yield by more than 95%, resulting in a negligible < 0.1% contamination. 

TABLE II: Systematic errors on determination of the branching fraction. 



Error Source Size (%) 

Finding the tt + track from the p + decay 0.3 

Hadron identification 1.0 

Finding the ir° from the p + decay 1.3 
-^extra < 0.2 GeV & 7T° efficiencies on background 1.1 

-E^extra < 0.2 GeV signal efficiency 2.0 

Background modeling 1.1 

Number of tags 2.0 

Tag bias 1.0 

Total 3^8 



Since we are requiring that -E ex tra be below either 0.1 or 0.2 GeV, it is necessary to 
check this efficiency in the data. Our procedure is to use the fully reconstructed sample of 
e + e _ — > D S D* S events selected the same way as described in Ref. |9[. The S ex tra distributions 
from Monte Carlo simulation and data for this sample are shown in Fig. [TUJ The agreement 
with the simulation is excellent. Table UTTI gives the efficiency from Monte Carlo simulation 
for our double tag sample and the efficiency measured in the data for specific ranges in 
-Eextra- The situation here corresponds to the extra energy deposited by two tags. We need 
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to translate these numbers to the case of one tag plus a signal p + . A cut at -E extr a = 0.3 
GeV in the double tag data corresponds to the same efficiency as a cut of 0.2 GeV in the 
single tag plus p + data. The difference between Monte Carlo simulation and data then is 
(—1.2 ± 1.6)%, which implies a systematic error of ±2.0% in this efficiency. 
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FIG. 10: The -Eextra distributions from e + e — ► D*D S events with both D s decays fully recon- 
structed (DT) for data (diamond) and Monte Carlo simulation (histogram). 



TABLE III: The efficiencies (e) for Data and Monte Carlo simulation for different requirements on 
E ex tra., and the corresponding fractional differences. 



-Eextra (GeV) £Pata(%) CMc(%) EPata/eMC ~ 1 (% 

<0.1 40.24 ± 1.27 40.81 ±0.31 -1.4 ±3.2 

<0.2 57.75 ± 1.28 59.12 ±0.31 -2.3 ± 2.2 

<0.3 72.35 ± 1.16 73.21 ±0.28 -1.2 ±1.6 

<0.4 83.27 ±0.97 82.91 ±0.24 0.4 ±1.2 



We note that if we fix the background branching fractions to their nominal values, and 
refit the data, the statistical error in the -E ex tra < 0.1 GeV bin decreases from 16.5 to 15.9 
events, and in the 0.1 < -E ex tra < 0.2 GeV bin decreases from 11.3 to 11.1 events. Thus, our 
statistical error contains a significant component from the background estimates. 
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B. Cross Checks Using 7r+7r° Helicity and Mass Distributions 



In principle the best way to view the p + polarization is to look at the angle 9 of the ir + 
with respect to the p + direction in the r + rest frame. Since we cannot reconstruct the r + , 
we use the laboratory frame. We consider all events in the -E ex tra interval below 0.2 GeV 
and having —0.05 <MM 2 < 0.60 GeV 2 . Fig. [TT1 shows cos# from the data with sidebands 
subtracted compared with the sum of expected signal and backgrounds from Monte Carlo 
simulation, normalized by yields from the data fit. The predicted total is in good agreement 
with the shape and data yield. 
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FIG. 11: Helicity distribution from the p + decay as measured in the laboratory frame. The points 
with error bars are the sideband subtracted data. The dashed line represents the predicted signal 
shape and the dotted line the predicted background shape from real decays. The Monte Carlo 
predictions are normalized by the fitted yields to the data. The predicted total is given by the solid 
line. We require that E extrSL < 0.2 GeV, and -0.05 <MM 2 < 0.60 GeV 2 . 

We also show the n + ir° mass distribution in Fig. [121 Again the predicted sum has good 
agreement with the shape and data yield. 



IV. CONCLUSIONS 

We list the CLEO-c measurements of leptonic branching ratios and fo s in Table IIVI To 
extract the decay constant we use M D + = 1.96849(34) GeV, a Df lifetime of 0.500(7) ps, 
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FIG. 12: Distribution of ir + n mass. The points with error bars are the sideband subtracted 
data. The dashed line represents the Monte Carlo predicted signal shape and the dotted line the 
Monte Carlo predicted background shape from real Df decays. The Monte Carlo predictions are 
normalized by the fitted yields to the data. The predicted total is given by the solid line. We 
require that E extla < 0.2 GeV, and -0.05 <MM 2 < 0.60 GeV 2 . 



and 1.77684(17) GeV for the r + mass |4|. While it has been customary to take |KJ = \V u d\ 
the expansion of the Wolfenstein parametrization of the CKM matrix to order A 4 |l3j implies 



that \Vcs\ = \V ud \ - \V cb \ 2 /2. We use \V ud \ = 0.97418(26) as derived in Ref. [ij. For \V cb \ we 
use a value of 0.04 from an average of exclusive and inclusive semileptonic B decay results 
as discussed in Ref. [15j . Thus, we find \V CS \ = 0.97338(26). The resulting value of fo a is 0.2 
MeV larger than taking \V CS \ = \V u d\. 

Previously reported values of fr> s have been corrected to corresponded to the above num- 
bers. These quantities contribute additional small amounts to the systematic error of ±1.8 
MeV (lifetime), ±0.1 MeV (V^), and for the t + v mode only ±0.4 MeV (M D +) and ±0.2 
MeV (t + mass), that are included in the quoted values. A theoretical upper bound on fq„ 



of 270 MeV has been calculated using two-point correlation functions by Khodjamirian [17 
The CLEO-c values for both the r + v and fi + u modes are below this limit. 

The ratio of decay constants from CLEO for the two leptonic decay modes is 

- — — — = 1.01 ±0.05, 6 

consistent with lepton universality. The average value of the pseudoscalar decay constant 
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TABLE IV: Recent absolute measurements of fjj g from CLEO-c 



Experiment 


Mode 


B (%) 


f Ds (MeV) 


This result 


T + V {p + V) 


(5.52 ±0.57 ±0.21) 


257.8 ± 13.3 ±5.2 


CLEO-c [9] 




(6.42 ±0.81 ±0.18) 


278.0 ±17.5 ±4.4 


CLEO-c [16] 




(5.30 ±0.47 ±0.22) 


252.6 ±11.2 ±5.6 


Average 


T + U 


(5.58 ±0.33 ±0.13) 


259.7 ±7.8 ±3.4 


CLEO-c [9] 




(0.565 ± 0.045 ± 0.017) 257.6 ± 10.3 ± 4.3 


Average 






259.0 ±6.2 ±3.0 



using both leptonic decay modes is 

f Ds = (259.0 ± 6.2 ± 3.0) MeV. (7) 

There are two SM based theoretical predictions for fu a in the literature based on Lattice 
QCD calculations, where all three light quark loops are included. The values predicted are 
(241 ± 3) MeV from the HPQCD±UKQCD collaboration [2], and (249 ±11) MeV from 
the FNAL±MILC±HPQCD collaboration [3]. We choose to compare with the more pre- 
cise Follana et al. result, realizing that it needs confirmation, especially with respect to 
the rather small error. The difference between the experimental average of fr> s and the 
HPQCD±UKQCD prediction is 2.4 standard deviations. Other theoretical predictions are 
given in Ref. [3j. 

Belle has also measured the absolute branching fraction for — > fi + u and found a 



value fo 3 = (275 ± 16 ± 12) MeV [19l . l20j ]. Combining with the CLEO measurements we 
find f Ds = (260.7 ± 6.5) MeV, which differs from the HPQCD±UKQCD prediction by 2.8 
standard deviations. We emphasize that this difference is qualitatively different than looking 
for new physics as a bump in mass spectrum where any values of the mass and width can 
be entertained. Here we are dealing with a theoretical number that was predicted before 
the measurements were available. Thus, although we cannot claim to have definitely seen 
an important discrepancy pointing to physics beyond the standard model, it is important to 
pay attention to this difference and to see what it may imply. 

In fact this possible discrepancy has motivated several new beyond the SM theories. 
These include leptoquark models of Dobrescu and Kronfeld [2lj|, R parity violating models 
of Akeroyd and Recksiegel [22|], and Kundu and Nandi who relate this discrepancy with 
preliminary indications of a large phase in 5c — B s mixing, and explain both with a specific 
supersymmetry model [23]. Dosner et al. [24] show however, that scalar leptoquark and R- 
parity violating models would have different effects on t + v and /i + z/ final states. Gninenko 
and Gorbunov argue that the neutrino in the D s decay mixes with a sterile neutrino, which 
enhances the rate and also explains the excess number of low energy electron like events in 



the MiniBooNE data 25 



We also have measured the following branching fractions: 



B(D+ -> K\ + tt°) = (1.00 ±0.18 ±0.04)% 



o, 



B(D+ -> tt+ttV) = (0.65 ±0.13 ±0.03)% 



o, 



B(DJ^rip + ) = (8.9 ± 0.6 ± 0.5 r... 
The first two modes have not been measured previously. 
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Appendix A: Monte Carlo Generated Background Modes 

We list the different modes that populate the MM 2 distribution for the three E extra inter- 
vals as given by the generic Monte Carlo simulation in Table EJ 

Appendix B: Measurements of Z)+ Branching Fractions for Selected Background 
Modes 

The K°tt + 7i° mode has not been previously measured. We select events opposite our D~ 
tag candidates with a single charged track consistent with being a tc + in conjunction with 
a 7T° candidate as described above, and a K s — ► 7r + 7r~ candidate where the invariant 7r + 7r _ 
mass is within 12 MeV of the known K s mass and the flight significance, the distance that 
the Kg travels divided by the error in the distance, is greater than 2. The invariant mass of 
i^7r + 7r° combinations is shown in Fig. [131 The data are fit with a signal CB function with 
all parameters except the area fixed to those given by the Monte Carlo simulation of this 
mode, and a second order Chebyshev background polynomial. The fit yields 44±8 events. 

We proceed by performing an unbinned liklihood fit to the Daltiz plot shown in Fig. [HJ 
using the isobar model formalism as described in Ref. [26]. The fit results are that the K s p + 
fraction is (88±8)% with only (17±8)% of A* + vr°; the relative phase is (21±25) degrees. The 
efficiency is determined to be 21% by Monte Carlo simulation which uses a Kp resonance 
structure, resulting in a branching fraction 

B(D+ -> A^tt+tt ) = (1.00 ± 0.18 ± 0.04)%, (9) 

where the systematic error arises from several sources shown in Table FVH The error due to the 
Dalitz plot structure is is 0.7%, found by evaluating the relative efficiency difference between 
pure Ksp + and the model resulting from our Dalitz plot fit; this is negligible compared to 
the other sources. Note that for our purposes some of the systematic error cancels because 
we are using the same D~ tag sample and the same 7r + and tt° detection efficiencies as for 
signal r + — ► p + z7. 

The 7r + 7T 7r° mode also has not previously been measured, though the analogous isospin 
related mode 7r + 7r + 7r~ has been. Here we require that the invariant 7r°7r° mass be more than 
50 MeV from the K s mass in order to reject K s . The invariant mass plot is shown in Fig. [TBI 
The signal CB function is fixed to the Monte Carlo predicted shape, and we use a second 
order Chebyshev polynomial function to model the background shape. We find a signal of 
72 ± 16 events. 
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TABLE V: Assumed branching fractions, Numbers of events (#) and fractions resulting from a 
generic Monte Carlo simulation of the Df backgrounds for twenty times the data in the interval 
- 0.2 < MM 2 < 1.0 GeV 2 . 

Mode B(%) < £ cxtra < 0.1 GeV 0.1 < £ cxtra < 0.2 GeV 0.8 GeV < £ cxtra 







# 


Fraction (%) 


# 


Fraction (%) 


# 


Fraction 


K°TT + TT° 


0.85 


348 


29.3 


128 


10.3 


548 


1.61 




7.58 


114 


9.6 


215 


17.4 


17797 


52.4 


TT+TT ^ 


0.58 


48 


4.1 


27 


2.2 


991 


2.9 


t + v, r + — ► 7r + 7r 7r°P 


0.55 


159 


13.4 


266 


21.5 


99 


0.3 


T + U, T + — > 7T + Z/ 


0.66 


63 


5.3 


18 


1.5 


2 


0.01 


t + u, t + — > other 


3.27 


81 


6.8 


43 


3.5 


12 


0.04 


V^" 3 V L b 


1.38 


70 


5.9 


114 


9.2 


560 


1.7 


X/j, + u 


5.87 


91 


7.7 


158 


12.8 


653 


1.9 




1.54 


15 


1.3 


32 


2.6 


1250 


3.7 




0.61 


19 


1.6 


15 


1.2 


9 


0.03 




3.67 


10 


0.84 


15 


1.2 


2610 


7.7 


K°K + ir° 


2.50 


20 


1.7 


7 


0.6 


32 


0.09 


K°K+ 


2.93 


11 


0.9 


6 


0.5 


31 


0.09 


K°7T + 

11 /I 


94 


q 


8 

VJ. o 


q 


7 


59 


1 7 




0.70 








1 


0.1 


156 


0.46 




1.25 


29 


2.5 


62 


5.0 


144 


0.42 




0.70 


18 


1.5 


8 


0.7 


31 


0.09 


Xe+v 


6.19 


8 


0.7 


5 


0.4 


16 


0.05 


K°7T+7T°TT 


0.65 


24 


2.0 


47 


3.8 


457 


1.34 


rjir + TT°TT 


3.25 


5 


0.4 


12 


1.0 


4357 


12.82 




3.87 


2 


0.2 


6 


0.5 


2041 


6.00 


a;7r + 


0.25 


1 


0.1 








18 


0.05 


7r+7r 7r°7rV 


0.85 


1 


0.1 


3 


0.2 


1684 


4.95 




7.35 


4 


0.3 


14 


1.1 


340 


1.00 


Other 




36 


3.0 


27 


2.2 


97 


0.29 



We also perform an unbinned liklihood fit to the Daltiz plot shown in Fig. [161 again using 
the isobar model formalism. There is no evidence for p + 7r°. We find that the fractions of 
/o(980)7r+, / 2 (1270)tt + and /o(1370)tt + are (56.5±9.1)%, (20.5±7.3)%, and (38.1±8.6)%, re- 
spectively. Fixing the / (980)7r + phase at zero degrees, the relative phases of the / 2 (1270)7r + 
and / (1370)7r + with respect to zero are (243±29) degrees and (299±24) degrees, respec- 
tively. 

The Monte Carlo simulated efficiency is 28.1%, yielding 

B(D+ -> tt+ttV) = (0.65 ± 0.13 ± 0.03)%. (10) 

The systematic errors are listed in Table [VTl This number is consistent with B(tt + tt + 7t~)/2 = 
(0.56±0.04)% 0, which is what is expected if the neutral dipion system final state dominates 
both modes. 
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FIG. 13: The invariant mass spectrum of Kgir + Tr . The curves show a second order Chebyshev 
polynomial function that describes the background summed with a signal CB function whose width 
is fixed (solid line). 

The r)p + branching fraction has been previously measured as (13.1 ± 2.2)% [4]; we wish 
to improve on this accuracy. We look for events with only one charged track consistent with 
being a ir + The 77 is looked for in the 77 decay mode only; mass combinations are used if 
they are within 3 standard deviations of the r] mass. We insist that the 7r + 7r° invariant mass 
be within 250 MeV of the p + mass. The resulting 777r + 7r° invariant mass spectrum is shown 
in Fig. [T71 Here we have enough data to let the r.m.s. width of the CB function vary in the 
fit. The background is again described by a second order Chebyshev function. We find a 
total of 328 ±22 events. We use a Monte Carlo determined efficiency of 22.4%. We find that 

B{Dj -> rjp + ) = (8.9 ± 0.6 ± 0.5)%. (11) 

The systematic errors are listed in Table [VTl Our new measurement is lower by about 1.8 
standard deviations than the PDG average [4|. 
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FIG. 14: (a) Dalitz plot of Kgir^ir , where the invariant mass of the three particles is selected 
within ±20 MeV of the Df mass, (b)-(d) show the mass projections, the solid curves show the 
overall fit and the dashed curve the background from sidebands. 



TABLE VI: Systematic errors on determination of the branching fractions of several background 
modes. We give two errors, one for the branching fraction to be used as an independent measure- 
ment (Ext), and the second is the one to be used internally (Int) for the lepton branching fraction 
analysis where some of the errors cancel. 



D+ K 7r+vr Df -» r]p + D+ -» vr+vrV 

Error Source (Ext)(%) (Int)(%) (Ext)(%) (Int)(%) (Ext)(%) (Int)(%) 



Hadron identification 


1.0 





1.0 





1.0 





Finding n° from {p + ) 


1.3 





1.3 





1.3 





Background modeling 


2.0 


2.0 


1.9 


1.9 


2.0 


2.0 


7T° efficiency 


1.3 





1.3 





2.6 


1.3 


K° efficiency 


2.0 


2.0 














rj efficiency 








4.0 


4.0 








Number of tags 


2.0 





2.0 





2.0 





Tag bias 


1.0 





1.0 





1.0 





Total 


3.8 


2.8 


5.1 


4.4 


4.0 


2.4 
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FIG. 15: The invariant mass spectrum of tt + tt tt candidates. The curves show a second order 
Chebyshev polynomial function that describes the background (dashed line) summed with signal 
CB function (solid line) whose width is fixed to the Monte Carlo predicted shape plus the measured 
Gaussian smearing. 
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FIG. 16: (a) Dalitz plot of 7r + 7r°7r°, where the invariant mass of the three particles is selected within 
±24 MeV of the mass. There are two entries per event, (b)-(d) show the mass projections, the 
solid curves show the overall fit and the shaded region the background from sidebands. The zero 
at 0.25 GeV 2 in (d) results from the K$ rejection criteria. 
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FIG. 17: The invariant mass spectrum of rj7r + ir candidates, where rj — > 77. The curves show 
a signal CB function whose width is allowed to float (dotted line), a second order Chebyshev 
polynomial function that describes the background (dashed line), and the sum (solid line). 



25 



[1] D. Silverman and H. Yao, Phys. Rev. D 38, 214 (1988). 

[2] E. Follana, C. T. H. Davies, G. P. Lepage, and J. Shigemitsu (HPQCD and UKQCD Collab- 
orations), Phys. Rev. Lett. 100, 062002 (2008). larXlv: 0706 .17261 [hep-lat]. 

[3] J. L. Rosner and S. Stone in Ref. [J| , see also larXiv:0802.1043l [hep-ex] . 

[4] C. Amsler et al. (Particle Data Group), Physics Letters B 667, 1 (2008). 

[5] A.G. Akeroyd and F. Mahmoudi, JHEP 0904, 121 (2009). larXTv:Q902.2393l [hep-ph]. 

[6] D. Peterson et al, Nucl. Instrum. and Methods A 478, 142 (2002); Y. Kubota et al. (CLEO 
Collaboration), Nucl. Instrum. and Methods A 320, 66 (1992). 

[7] M. Artuso et al, Nucl. Instrum. Methods A 554, 147 (2005) [physics/0506132| . 

[8] D. Cronin-Hennessy et al. (CLEO Collaboration), Phys. Rev. D 80, 072001, (2009), 
larXiv:0801.3418l [hep-ex]. 

[9] J. P. Alexander et al (CLEO Collaboration), Phys. Rev. D 79, 052001 (2009). la7XTv:0901.1216l 
[hep-ex] . 

[10] M. Fujikawa et al (Belle Collaboration), Phys. Rev. D 78, 072006 (2008), arXiv:0805.3773l 
[hep-ex] . 

[11] G. Bonvicini et al. (CLEO Collaboration), Phys. Rev. D 70, 112004 (2004), 

|arXiv:hep-ex/0411050~] 
[12] T. Skwarnicki, DESY F31-86-02 (thesis, unpublished) (1986). 
[13] J. Charles et al. Eur. Phys. J. C 41, 1 (2005), [arXTv:hep-ph/0406184 
[14] I. S. Towner and J. C. Hardy, Phys. Rev. C 77, 025501 C2008) larXw:0710.3181l [nucl-th]. 
[15] M. Artuso, E. Barberio and S. Stone, PMC Physics A, 3:3 (2009) larXw:0902.3743l [hep-ph]. 
[16] P.U.E. Onyisi et al. (CLEO Collaboration), Phys. Rev. D 79, 052002 (2009). la7Xiv:0901.1147l 

[hep-ex] . 

[17] A. Khodjamirian, Phys. Rev. D 79, 031503(R) (2009) la7Xiv:0812.3747l [hep-ph], 

[18] C. Bernard et al. (Fermilab Lattice, MILC and HPQCD Collaborations), arXiv:904.1895 

[hep-lat] (2009); C. Aubin et al, Phys. Rev. Lett. 95, 122002 (2005), [a7xTv:hep-lat/0506030 
[19] L. Widhalm et al. (Belle Collaboration), Phys. Rev. Lett. 100, 241801 (2008). larXTv:Q709.134QI 

[hep-ex] . 

[20] We have not included other measurements of fu g in the average that are based on measuring the 
leptonic branching fractions with respect to a hadronic decay mode due to the relatively large 
uncertainities in the determination of the relevant hadronic branching fraction. See Ref. 0] for 
a compilation of these results. 

[21] B. A. Dobrescu and A. S. Kronfeld, Phys. Rev. Lett. 100, 241802 (2008), arXiv:0803.0512l 
[hep-ph] . 

[22] An example of such a model is R-parity violating SUSY; see A. G. Akeroyd and S. Recksiegel, 

Phys. Lett. B. 554, 38 (2003), [hep^ph/0210376[ 
[23] A. Kundu and S. Nandi, Phys. Rev. D 78, 015009 (2008). larXTv:Q803.1898l [hep-ph], 
[24] I. Dorsner, S. Fajfer, J. F. Kamenik, and N. Kosnik, arXiv:0906.5585 [hep-ph]. 
[25] S. N. Gninenko and D. S. Gorbunov. larXiv:0907.4666l [hep-ph] (2009). 

[26] S. Kopp et al. (CLEO Collaboration), Phys. Rev. D 63, 092001 (2001); P. Rubin etal (CLEO 
Collaboration), Phys. Rev. D 78, 072003 (2008). 



26 



